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INTERIOR  BALLISTICS  OF  RECOILLESS  GUNS 


Abstract 

A  system  of  interior  ballistics  is  developed  for  guns 
in  which  the  recoil  is  eliminated  by  the  rocket  action  of 
powder  gas  flowing  through  a  Venturi  in  the  breech.  The 
system  is  illustrated  and  applied  to  numerous  examples  in¬ 
cluding  the  German  7-5-cm  Leichtes  Geschtitz. 

1 .  Introduction 

In  this  report  we  present  a  study  of  the  interior  ballistics  of 
recoilless  guns.  The  equations  for  the  interval  of  burning  of  the 
powder  are  reduced  with  the  help  of  simple  approximations  to  the  equa¬ 
tions  valid  for  the  corresponding  interval  in  conventional  guns.  An 
approximate  solution  for  the  interval  after  the  powder  has  burned  is 
readily  obtained  when  it  is  noticed  that  the  adiabatic  law  holds  in 
this  case  as  well  as  in  conventional  guns  and  rockets.  The  accuracy 
of  the  approximations  is  checked  by  comparing  a  solution  with  the  re¬ 
sult  of  a  numerical  integration  of  the  fundamental  equations. 

The  ballistic  system  developed  for  recoilless  guns  makes  possible 
prediction  of  the  performance  of  such  weapons  with  an  accuracy  that 
should  be  comparable  with  the  results  obtained  by  any  of  the  leading 
ballistic  systems  for  conventional  guns.  The  methods  presented  here 
can  bo  used  to  estimate  the  performance  of  any  gun  with  a  gas  leakage 
through  a  vent  of  constant  area,  such  as  smoothbore  mortars  or  vented 
mortars  (where  the  muzzle  velocity  is  controlled  by  varying  the  gas 
leakage.)  In  other  cases  of  gas  leakage  where  the  vent  area  is  not 
constant  throughout  the  travel  of  the  projectile  —  such  as  worn  guns 
or  guns  with  gas  operated  automatic  devices  —  it  is  necessary  to 
integrate  our  fundamental  equations  numerically. 


PART  I.  THE  RECOILLESS  GUN 


The  recoilless  gun  is,  in  essence,  a  gun  with  an  abnormally  large 
chamber  that  is  terminated  in  the  rear  by  an  open  nozzle  or  Venturi  in 
place  of  the  conventional  breech  block  (Fig. l).  Thus  the  recoiling 
forces  are  balanced  by  the  rocket  action  of  powder  gases  escaping  to 
the  rear.  The  Germans  are  using  this  type  of  gun  at  the  present  time. 


Fig.  1.  Schematic  diagram  of  a  recoilless  gun. 

and  the  English  are  experimenting  -with  several  models.  We  believe  that 
recoilless  guns  may  well  be  .superior  to  rockets  and  to  conventional-  guns 
for  a  number  of  tactical  uses.  The  main  disadvantage  of  such  guns  is, 
of  course,  the  blast  of  gases  t,o  the  rear.  This  restricts  their  use  to 
tactical  situations  in  which  a  danger  area  to  the  rear  of  the  gun  may  be 
kept  free. 

3 .  The  lack  of  recoil 

Fortunately,  the  suppression  of  recoil  does  not  require  a  delicate 
balancing  of  forces.  In  fact,  if  the  area  of  the  throat  of  the 
nozzle  is  the  same  as  the  cross-sectional  area  A  of  the  boro  and  if  the 
nozzle  has  no  flare  to  the  rear,  then  (neglecting  the  friction  of  the 
projectile  and  pressure  gradients  in  the  chamber)  there  is  no  resultant 
force  exerted  on  the  gun  by  the  powder  gasj  hence  there  can  bo  no  recoil. 
Actually,  because  of  the  friction  of  the  projectile,  there  would  be  a 
small  force  forward  when  A  and  A^.  are  equal.  In  addition,  the  gas 
escaping  at  high  velocities  through  the  nozzle  creates  a  pressure  drop 
along  the  back  wall  of  the  chamber,  so  that  there  is  a  further  forward 
thrust  exerted  on  the  gun  even  if  A  =  A^ .  Thus,  in  order  to  have  the 
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gun  completely  recoilless,  it  is  necessary  to  design  it  so  that  <  A. 
The  exact  size  of  A^  -will  depend  on  the  design  of  the  nozzle,  and,  in 
particular,  on  the  escape  velocities  reached  by  the  powder  gases.  If 
the  nozzle  is  so  designed  that  the  powder  gases  have  a  maximum  velocity 
of  the  order  of  7000  ft/sec,  then  a  value  of  approximately  1 . US  for 
A/A^  would  be  effective  in  eliminating  the  recoil.  In  the  German  75-mm 
rocoilless  gun,  the  ratio  A/A^  is  I.J48.  High  escape  velocities  lead 
to  an  efficient  gun  —  that  is,  a  small  powder  charge  —  but  high  escape 
velocities  also  increase  the  danger  zone  to  the  rear  of  the  gun. 


PART  II.  THE  BALLISTIC  EQUATIONS 
I4 .  Gas  flow  through  the  nozzle 

Let  N  (lb)  bo  the  weight  of  powder  burned  at  any  time  and  N'  (lb) 
the  weight  of  powder  gas  remaining  in  the  gun  at  any  time,  so  that 
d(N-N')/dt  is  the •  time-rate  of  flow  of  the  mass  of  gas  through  the 
nozzle.  Let  Tc  (°K)  be  the  isochoric-^  flame  temperature  of  the  powder 
and  T(°K)  the  actual  temperature  of  the  powder  gas  in  the  chamber. 
Finally,  let  P (lb/in?)  be  the  average  pressure  in  the  gun  at  any  time. 
Then,  according  to  the  usual  engineering  equation  for  gas  flow  through 
a  nozzle,-^ 

d(N  -  N')/dt  =  kVf”7TAtP  lb/ sec .  (l) 


The  nozzle  coefficient  k  may  be  found  from  the  ratio  of  specific 
heats’?  and  the  “impetus11  of  the  powder,  F  =  nRTQ  (ft  lb/Lb ) ,  where  n 
is  the  number  of  moles  of  gas  formed  from  unit  weight  of  powder  and  R 
is  the  gas  constant  per  mole.  The  equation  for  the  nozzle  coefficient 
k  is 


k 


/32.17U? 

/  F 


2  \(?  +i)/{?-i) 


sec-1 . 


(2) 


]_/  Here  the  word  isochoric  is  used  to  indicate  the  flame  tempera¬ 
ture  at  constant  volume. - 

2/  See,  for  example,  Stodola  and  Lowenstein,  Steam  and  gas  tur¬ 
bines  (McGraw-Hill,  1927),  vol.  I,  p.  '4I4,  Eq.  (15>). 
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For  FNH-M2  ponder  (Hercules,  20  percent  N.G.), 

k  =  0.005980  sec-1  [?  =  1.223  and  F  «  38UOOO  ft]. 

For  FNH-M1  powder  (du  Pont,  single-base), 

k  «  0.006718  sec”1  [>  =  1.258  and  F  =  310000  ft]. 

i  ... 

The  constants  for  any  powder  may  be  determined  by  the  methods  given  in 
NDRC  Report  A-li|2.^/ 


5.  Rate  of  burning  of  the  powder 

We  suppose  that  the  powder  burns  linearly  at  a  rate  which  is  pro¬ 
portional  to  the  pressure  and  has  a  constant  burning  surface.  Then 

H  =  C  |  P  Ib/sec,  (3) 

where  C  (lb)  is  the  total  weight  of  the  powder  charge,  W(in. )  is  the  web 
thickness  and  B (in./sec)/(lb/in?)  is  the  so-called  burning  constant .  If 
the.  powder  does  not  have  a  constant  burning  surface  it  is  necessary  to 
introduce  a  form  function  at  this  point. 


6 .  Rate  of  increase  of  powder  gas  in  gun 

Combining  Eqs.  (l)  and  (3),  wo  obtain  the  rate  of  increase  of  pow¬ 
der  gas. 


dN* 

dt 


cb 

w 

L- 


P  lb/sec. 


(2) 


We  assume  that  the  powder  gas  is  retained  in  the  chamber  by  a  blow-out 
disk  until  the  projectile  starts  to  move.  Let  NQ  be  the  weight  of  pow¬ 
der  burned  up  to  this  time.  It  is  convenient  to  define  9  by  the  rela¬ 
tion, 

_  f  -  k/pAt  _  _ 

CB/W  CB/V  ' 


3/  J.  0.  Hirschf elder,  R.  B.  Kershner  and  C.  R.  Curtiss,  Interior 
ballistics,  I,  NDRC  Report ,A-1U2  (0SRD  No.  1236).  This  report  will  be 
referred  to  hereafter  as  A— 1 U2. 
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It  mil  be  noticed  that  9  is  a  linear  function  of  \/l  /T.  In  terms  of  9, 


Eq.  (U)  may  be  written, 


9  ^  P  =  9  lb/sec 
W  dt  ' 


or,  after  integrating, 


where 


N'  =  9  N  +  (1  -  9  )N  lb, 
avg  v  avg'  o  ’ 

9a  is  an  average  value  of  9  during  the  burning. 

d  Vg  — 


7 .  Equation  of  motion 

The  average  pressure  P  in  the  gun  is  somewhat  larger  than  the  pres¬ 
sure  Pj  on  the  base  of  the  projectile  owing  to  the  existence  of  a  pres¬ 
sure  gradient  caused  by  the  inertia  of  the  powder  gas.  In  a  conventional 
gun,-/ 

P  -  Win?,  (6) 

M 

where  M  (lb)  is  the  weight  of  the  projectile  and  S  is  a  parameter  depend¬ 
ing  on  C/k  and  having  a  value  a  little  larger  than  3  in  most  cases.  In 
a  recoilless  gun,  since  the  over-all  mcmentums  are  balanced,  there  is  a 
position  of  zero  gas  velocity  in  the  chamber  and  it  seems  reasonable  to 
assume  that  the  pressure  gradient  established  forward  of  this  position 
of  zero  velocity  is  about  the  same  as  that  in  a  conventional  gun.  Thus 
we  assume  that  a  relation  similar  to  Eq.  (8)  is  valid,  the  single  dif¬ 
ference  being  that  only  the  fraction  of  the  charge  that  is  retained  in 
the  gun  should  be  used.  If  NQ  is  small,  then  it  is  clear  from  Eq.  (7) 
that  the  charge  weight  retained  in  the  gun  is  approximately  Q^vp^* 

Thus  we  assume  for  a  recoilless  gun, 

P  =  FX(M  +  9aycrC/S)/M  lb/in?  (?) 


It  Is  convenient  to  express  the  equation  of  motion  in  terms  of:  the 


h/  Report  A-lli.2,  Eq.  (2). 
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average  pressure  P(lb/in?)j  the  velocity, 

V  ’  T2  I  ft/sec> 


(10) 


where  X(in. )  is  the  effective  distance  frcm  breech  to  projectile 
defined  so  that  AX  (in?)  is  the  total  bore  volume  behind  the  projec¬ 
tile  ?  and  the  effective  mass  of  the  projectile, 


m«  -  1.0U  (M  +  eavgC/3)/32.17  slugs 


or,  to  a.  sufficiently  good  approximation, 


m'  =  I.OI4 


M  +  1  - 


kA. 


/ 32.17,  '• 


(11) 


(IT) 


C(B /V)J  31 

where  the  factor  1.0U  accounts  for  friction.  In  terms  of  these  quan¬ 
tities,  the  equation  of  motion  becomes 


p={ILLfiZ=i2~  V—  lb /in2 
A  dt  1  A  dX 


(12) 


8.  The  temperature  in  the  chamber 


•The  internal  energy  of  the  gas  in  the  gun  is  given  by  the  equa¬ 


tion, 


N* 


f  C  dT  =  N  / 
J0  —  ^ 


/’N-N1  /.T 


cv  dT  -  / 

0  —  "0  Jq 


^dT  d(N-N')  -Im’V2  -/slm'V?  (13) 


The  notations Cy  and  indicate  the  specific  heat  at  constant  volume 
and  constant  pressure,  respectively,  as  functions  of  the  temperature. 
Here  N 1° 


/•N-N'  ? 

Jo  Jo 


-q  Gv  dT  is  the  total  energy  of  the  powder  gas  that  is  formed? 
T  — 


/0  dT  d(N-N')  is  the  kinetic  and  thermal  energy  of  the  gas 

that  has  escaped  [it  will  be  remembered  that  1  lb  of  powder  gas  ex¬ 
panding  adiabatically  through  a  nozzle  into  a  region  of  low  pressure 

acquires  the  kinetic  energy,  -|V2/32.17  -  C  dT  ft  lb]?  -jgm'V2  is  the 

</0  — i- 

kinetic  energy  of  the  projectile  and  of  the  powder  gas  in  the  gun,  to¬ 
gether  with  the  work  going  into  friction?  and  P>\ m'V2  is  the  heat  loss 
frcm  the  powder  gas  to  the  bore  (assumed  proportional  to  m'V2). 
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If  we  set  C!  =  C,r  +  nR,  let  /  °C_r  dl  =  En  and  let  C  be  the  aver- 
_P  _V  (/q  _Z  o  v 

age  value  of  In  the  high  temperature  range  between  T0  and  T,  then 
Ecf,  (13)  becomes 


rN-N  * 


N'[E  -C  (T  -T)]  =  NE  -  [nRT+S  -C  (T  -T)]  d(N-N')-  (l  +  ^)|m'V2 

uov  o,  o  J  n  ovo  ' 


»  NE 


/>N-N 1 

-  (Eq- CvTq)-(N- N ')- J  (nR  +  Cy)T  d(N-  N«)-  (1  +/j)|m'V?  (lU) 


Collecting  terms  and  setting  }  =  (riR  +  C  )/C  ,  we  obtain 

V  V 


N»C  T  =  NC  T  -  ?C  , 

V  VO  V  J 


[ 


N-N» 


0 


T  d(N  —  N* )  -  (1  +  (i)2i  m‘V? 


OS) 


If  -g-m'V2  =  0  and  NQ  =  0  (so  that  N'  =  ©N)  the  steady  temperature 
solution  of  Eq.  (IS)  is  given  by  the  equation. 


To  =  "  9('}  ~  l)]*  06) 

Thus,  if  9  =  1 ,  corresponding  to  a  closed  chamber,  then  =  Tq,  the 
isochoric -flame  temperature.  Of  if  ©  =  0,  corresponding  to  the  steady 
state  in  a  rocket  motor,  then  =  T  /?,  the  isobaric  flame  temperature. 

In  order  to  calculate  the  energy  of  the  powder  gas  before  all  of  . 
the  powder  is  burned,  it  is  necessary  to  approximate  the  integral  cor¬ 
responding  to  the  energy  of  the  gas  that  is  discharged.  Making  use  of 
Eqs.  (l)  and  (12),  we  have 


d(N-  N1 ) 
dt 


dt=  >C 

v 


k  A  i/f 


oxavg 


!*  t 


Pdt  =  0-C  k  \ZTt 

v  0  avg 


(A  /A)m‘V. 


(17) 


Here  T  is  a  suitable  average  temperature  in  the  chamber  during  the 
gas  discharge.  This  integral  has  been  evaluated  numerically  for  hypo¬ 
thetical  recoilless  guns  and  0?  obtained  as  a  function  of  the  veloc- 
_ avg 

ity.  At  first  vF  equals  vF~,  then  drops  rapidly  and  thereafter  is 

c,  vg  o 

linear  with  velocity.  We  make  very  little  error  if  we  set 


\/f 

avg 


er 


n/F 

0 


id-^7W)^] 


(18) 


-  8  - 


where 


Ttt  =  T  / 
o  o/ 


kA  - 

1  +  — _t —  ()  -  1 ) 

C(B/W) 


(19) 


Here  Tb  and  Vb  are  the  temperature  and  velocity  at  the  time  when  all  of 
wder 


the  powder  is  burned.  With  this  definition  of  yT_  ,  Eq.  (17)  becomes 

ca.  v  0 


v 


'0 


T  d(N  -  N< )  -  ?Cvk  (-fJm'V  -  |(1  -  ^A7^) 


V 


.  (20) 


If  wc  define  }  as  in  A— 1 U2 ,  by  the  relation, 

T  -  1  *  (1  +  ,£)(>  -  1), 


(21) 


and  let 


cvto  =  F/(y  *r  1). 


(22) 


the  energy  equation  is  expressed  in  terms  of  T  and  V  by  the  relation, 
T 

•  =  N - - - - - 7  k 

A  if  T 


N» 


(>-  uin'v2  «t  K  r  ■•(  ny\  v 

0  o  b 


.  (23) 


9 •  The  equation  of  state 

Let  A  (lb/in?) ,  the  density  of  the  powder  gas  at  any  time,  be  given 
by  the  equation,  ... 


N» 

•  AX  -  (C  -  N)//0  * 


(-2U) 


Here  AX  is  the  volume  behind  the  projectile j  and  p  is  the  density  of  the 
solid  powder  so  that  (C  -  N )/p  is  the  volume  occupied  by  solid  powder. 
For  the  equation  of  state  of  the  powder  gas,  we  use  the  Abel  form,  with 
the  covolume  v^(in?/lb) , 


12F 


I'* 


P_ 

N' 


AX  -  ~J!  -  7N']  in. 


J 


(25) 


(in?),  the  volume  of  the  chamber,  and  A  =  C/vc  (lb/in?) , 


Now  AXq  15  V  g  \  J-ii*  J  y  UiiC  VVJ.UIU^  WJ.  uxiv?  1  If^XU  J.  y  wnw.  ^ Q  w/  VQ 

the  initial  density  of  loading.  Thus  Eq.  (25)  becomes 


N’T®  -S-  v 
12F  c 


X 

Lxl 


P 


A_  N 

P  c 


+  -  -  r?A  —  P  lb°K. 

'  O  r> 


N» 

~cj 


(26) 
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In  the  small  covolume'  correction  term  7^-g-we  approximate  N'  by  the 

use  of  the  constant  9  from  Eq.  C 7)  S  that  is, 

avg 

N»  =  9  N  +(1  -  9  )N  .  (2r 

avg  avg  o 

Also  N  may  be  found  by  integrating  Eq.  (3)  "with  the  help  of  Eq,  (12); 
this  gives 

C  G  AW 

Substitution  of  Eqs .  (27)  and  (28)  into  Eq.  (26)  and  introduction  of 

a  =  y[  -  (1  / f) 


V  -  (!//>) 
3s  i  -  040 


gives 


T  f  X  A  N  ,  m  ’  B  *[ 

N»  T  =  -2-  v - 2-aA-~-a  A  ™  v 

1?F  CY  n  OQ  0  AWJ 


12F  C|X 


Here  j*  plays  the  same  role  as  j3  in  A-1U2. 


1 0 .  The  ballistic  equation; 


According  to  Eq.  (28), 


N  W0  m1  B 

-  -  -£  + - V. 

C;  .  C  ,A.W 


In  particular,  if  denotes  the  velocity  at  the  instant  when  the  pow¬ 
der  is  completely  burned,  that  is,  when  N  =  C,  we  have 

v  «  1  “.VP—..  (33) 

b  (m'A)(BA0 

If  we  eliminate  N'T  from  Eq.  (23)  by  the  use  of  Eq.  (31)  and  N  by  the  use 
of  Eq.  (32),  there  results 

IS-  1  )m'v3  +  -  j  -  aAo  F  -  a4o33  r  s  v]p 

=  or  [^  +  rfv->ki5m'^v 
+  a  -v5v7t7)  ?]• 


AC  o'  o 


b'  o'  Vv 


-  10.  - 


Now  eliminating  P  by  the  equation  of  motion,  Eq,  (12),  and  using  Eq. (33) , 
we  get  the  fundamental  ballistic  equation  for  the  burning  interval  in 
the  form  of  a  differential  equation  for  the  velocity-travel  curve, 
namely, 


m* 

rx. 

Ao 

No 

f  m’ 

B  ‘ 

dV 

V2  + 

— 

V. 

-2  -  aAn 

—  aAJ. 

t  _ 

-  V 

V 

— 

A 

c 

k>  ’ 

P 

C 

0 

3  A 

W  J 

dX 

"N0 

B  ' 

An 1 

B\2  ~ 

+  • 

51  — 

-  V  +  ki 

— 

-  v2 

y 

=  CF^ 

.C 

1  A 

¥  3' 

\A 

¥/  J 

where 


and 


„  ,  ,  .  IkW'o 

C  ( B/lfj 


31 


k'  - 


jkA/X  (1  -  v^b/T^)' 


(35) 


(36) 


(37) 


"2  2 G(B/W)  (1  -  No/0) 

This  equation  is  seen  to  have  the  same  form  as  Eq,  (26)  of  A-1I42.  Fol¬ 


lowing  the  procedure  used  there,  we  introduce 

Ao  A  N0 
*  =  —  +  aAo  c~> 

_  GFm1  (B/¥)3 
61  A2 


_CF  B  .1 
e2  A  W 


u  •  §(*  -  i  )  -  k^  e-,_, 
N0/C 


q  = 


•  1  \  a  1 


ei(ji.)‘ 

r  =  aA0j^  j;ei, 

Z  «  V/ea,  ' 


(38) 

(39) 

(1*0) 

Oil) 

0l2) 

(1*3) 

0*1*) 


X  " 

y  “xt-01' 


0*5) 


and 


These  substitutions  reduce  Eq.  (35)  to  the  form, 

_  Z(y  -  rZ) 
dZ  q  +  Z  -  uZ2* 

which  is  identical  with  Eq.  (38)  of  A— 1 14.2 . 


(U6) 


Since  the  initial  conditions  used  here  are  the  same  as  in  A-Il(2, 
namely,  Z  =  0  when  X  =  X  ,  the  velocity- travel  relation  will  have  the 
same  form  as  Eq.  (1(8),  of  A-1i(2  —  that  is, 


(U7) 

(U8) 

(1*9) 

(50) 

(51) 


1 1 .  The  calculation  of  and  9pVg 


Tho  temperature  at  the  time  when  the  powder  is  completely  burned 
can  be  found  by  placing  V  =  V^,  N  =  C,  N1  =  and  T  =  T^  in  Eq.  (23). 
This  gives 


Tb  m  (F-D 


I V2  -  Ik—  v^Vf-  m '  V,  -  -  k  — 


m'VbV^*  «2) 


All  the  coefficients  in  this  equation  for  Vu^/Tq  are  easy  to  evaluate, 
with  the  exception  of  N^.  The  factor  cannot  be  found  easily,  and 
we  resort  to  an  approximation  method.  We  start  with  Eq.  (I4),  which  by 


Eq.  (12)  may  be  written, 

d??* 

dt 


-  12  - 


(53) 


’  B 
C  -  -  kl—  A*. 

.  W  V  T  ^ 


m*  dV 

j  T  dt ' 


We  assume  that  V^/t  is  a  linear  function  of  the  velocity  running  from 

V*T  /Tw  to  ]/r  7tT  so  that 
o'  o  o'  b 


V 


\/T  /T  =  /T"  ~  (v^T  /T!l  -  7T) 

o'  o  a  o  o  o  b  V, 


(5b) 


Then,  substituting  this  expression  into  Eq.  (53)  and  integrating,  we 
have 


Nf  Nq  m1  B 


kA+ 


m 1  B 


- {\fT~pF  +  \f5TFF) - v,  (55) 

C  C  A  W  b  2  C(B/W)  0  0  0  b  A  w  b 


-  V,.  - 


or,  making  use  of  Eq.  (33), 


N*  1  /  N  \  kA,  ,  ■  ,  ll’  r—r— 

f  -  1  -  i(i  -  -*)  (vryw  *  vC7tD. 


C  /  C(B/W) 


o'  o 


o'  b ' 


(56) 


Since,  according  to  Eq.  (5), 

9  *  1 

it  is  convenient  to  introduce 


k  V*Tq/T  At ' 

c(b/w)  * 


9”  =  1 


ky^n/T»  A+ 


C(B/W) 


Then  Eq. (55)  becomes 

K 


=  1  -  i(l  -  ^Jd  -  ©J)0  +  ^7*b)- 


(57) 


(58) 


(59) 


Substitution  of  Eqs.  (59)  and  (58)  and  (33)  into  Eq.  (52)  gives 


1-1(1- 


L  2 


(1-y)(i-9o)v5X 


Vf, /t  1 

b'  o 


1  -  I(1  -  ^)(1  -e")  I*-  w-”'7?. 


(60) 


The  temperature  when  the  powder  is  completely  burned  is  found  by  solv¬ 


ing  Eq.  (60),  a  quadratic  equation  in  VT^/T^ 


13 


It  is  found  that  a  very  satisfactory  value  of  ©avg  for  uso  in  con¬ 
nection  with  Eq. (27)  is  found  by  taking 


kv^V^avg  ^t 


9  '  -■■■ara-  , 

avg  c(3  fir) 


(61) 


where 


T  =  i(Tu  +  T,  ). 
avg  A  0  b 


(62) 


12,  The  functions  J  and  S 

The  velocity-travel  relation  involves  the  two  functions  J  and  S, 
given  by 

rZ 


J  ~  exp 


q  +  Z  -  uZ' 


dZ 


(63) 


and 


S  =  J  f  (1/J)  dZ. 
rf0 


(61t) 


Tables  of  these  two  functions  covering  the  range  that  usually  occurs 
in  ballistics  problems  are  included  in  A— 1 14.2 »  However ,  the  values 

of  q,  u  and  Z  that  occur  in  a  typical  problem  of  recoilless  gun  bal¬ 
listics  are  outside  the  range  of  the  tables  of  A-llj.2.  Thus  it  is 
desirable  to  have  a  convenient  method  for  computing  J  and  S.  Since 
the  integral  defining  J  is  elementary,  the  function  _J  can  be  calcul¬ 
ated  accurately.  However,  the  integral  defining  S  is  not  elementary 
so  it  is  necessary  to  resort  to  an  approximation  formula.  Fortunately, 
an  error  of  5  percent  in  S  is  permissible,  since  S  enters  the  ballistic 
equations  multiplied  by  a  coefficient  which  is  small  compared  -with  that 
of  J.  Let 


rx  =  (\/1  +  Uuq  -  l)/(2q). 

(65) 

r2  =  (\/l  +  Iraq  +  l)/(2q). 

(66) 

ti  “  l/(r1Vl  +  iinq), 

(67) 

t2  -  1  / (rs  Vi  +  iiuq). 

(68) 

-  Tli  - 


Then  a  convenient  form  of 'the  exact  formula  for  J  is  •' 

J  =  (1  -  r  1Z)~151  (1  +  r2Z)“ta.  (69) 

An  approximation  formula  for  S  usually  valid  to  -within  2  percent  is 

S=[j+(l-r1Z)  tl-(t-r1Z)[l +(l+raZ)“ta3}/[2n(t1 +1)].  (70) 


13.  After  the  powder  is  burned 


(a)  The  temperature  in  the  chamber.  —  After  all  of  the  powder 
is  burned,  the  weight  N*  of  gas  in  the  gun  becomes  a  decreasing  func¬ 
tion  of  time  while  N[=  C]  is  constant.  Let  and  T^  be  the  values 
of  N'  and  T  at  the  instant  the  powder  is  all  burned.  Then  the  equa¬ 
tion  of  energy,  Eq.  (1£>),  may  be  written  as 


N'C  T,  -  N'C  T  =  A 
b  v  b  v 


(71) 


It  will  be  shown  that  this  energy  equation  leads  to  the  adiabatic  ex¬ 
pansion  law  in  the  case  of  an  ideal  gas.  To  show  this  we  assume  the 
equation  of  state  of  an  ideal  gas,  that  is, 


PAX  =  N'nRT  in.  lb,  (72) 

and  the  adiabatic  law  connecting  the  pressure  and  density,  that  is. 


P(AX/N'/  «  K, 

where  the  constant  K  may  be  'evaluated,  frcm 


(.73) 


K  -  Pb(AXb/N^)’- 


(7U) 


if  is  the  pressure  at  the  instant  the  powder  is  all  burned  and  X^ 
is  the  value  of  X  at  the  same  instant.  From  Eqs.  (72)  and  (73)  the 
temperature  is  found  as 


-  JL  (EL)7"1 

"  nR  \  AX./ 


(75) 


It  will  now  be  shown  that  Eqs.  (73)  and.  (75)  provide  a  solution  of  the 
energy  equation,  Eq.  (71).  According  to  Eq.  (75),  the  last  term  of 


-  i5  - 


Eq.  (71 )'  becomes 


/>N '  c  K  /*N«  , 

/  C  T  dN'  «  -2-  /  (AX)~(* -1)  d(N'  )?  . 


4'  p 


J  nR  5/  AT 1 


Integrating  by-  parts  and  using  =  20^,  we  obtain  from  Eq.  (76), 


T  I  P 


C  T  dN'  =  (AX)“U-l)  (N»f  -  —  (AX.  (N,*  )? 

"  nR  D  n 


+  0-1)  A -2—  /  (N')’  (AX)“?dX.  .  (77) 

nR 

Substituting  this  value  of  the  integral  term  in  Eq.  (71)  and  elimina¬ 
ting  the  temperature  by  Eq.  (75) ,  we  get 


K  /N- 


0  v  nR  \  AX- 


1  \?-i 


-  N'C 


K  /N! 


nR  UX 


-A  /  PdX 


+  (N'f  +  (AXb)~(?_1^  (M^)7 


-  (7  -  1)  A 


(IV)7  (AX)“?dX, 


PdX=  0  ~  0 


AC„K  fX/N'O 


But  Eq.  (79)  is  satisfied  if  the  pressure  is  given  by  Eq.  (73) ,  since 

(1-1)0  “  nR. 

v 

(b)  Ballistic  equations.  —  The  equation  for  gas  flow  gives 


dN'/dt  =  -  d  (C  -  N*  )/dt  =  -k  Aty^7T  P. 


Fran  the  equation  of  motion, 


m '_av 

A  dt 


12m'  „  dV 
A  V  dX‘ 


Finally,  from  Eq.  (73), 


P(AX/N')?  =  K  = 


Pb(AXb/up» 


(82) 
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If  we  assume  that  the  temperature  T  occurring  in  Eq.  (80)  is  constant, 
then  the  ballistic  equations  for  the  interval  after  the  powder  is 
burned  may  be  obtained  from  Eqs.  (80),  (81 )  and  (82). 


In.  Eq.  (80),  let  \fiQ/ T  =  V^/T^  =  const.  Then,  from  Eqs.  (80)  and 


(81), 


dN«  A,  _____  dV 

Jk  ,  rv  /m  ~  r  __ 


dt 


m' 


dt 


or,  integrating, 


km'\/ro/  Th  /  \ 

N'  =  Nb  -  (v  "  V 


(83) 


(81*) 


Let 


MJ,  (&/At) 


and 


v  =  V/Vh. 


Then  Eq.  (81|)  may  be  written, 


N*  =  N£[1  +  $  -  jrfv]. 


(83) 

(86) 

(87) 


Now  P  and  N'  may.be  eliminated  from  Eq.  (82)  by  the  use  of  Eqs.  (81 ) 
and  (87).  This  gives 


12m'  dV 

.  _  T7 

AX 

7 

=  P 

I 

IcT 

_ j 

A  dX 

uN^(l  +  4  -  ^v)_ 

*b 

ln;  J 

Using  Eq.  (86)  and  introducing 

1  -  ■  - 

-APl.  V’/IZ" ''Vb  « 

and 


7  * n 


we  may  write  Eq.  (88)  as 


dv 

dY  V  V 


Cl. 


(88) 

(89) 

(90) 

(91) 

(92) 
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Eq.  (92)  may  be  integrated  by  separating  variables.  The  result  is 


SI 


dY 


r 


Ji 


vdv 


(V  -  v) 


(93) 


or 


11(2  -  ?)|  1  -  — : 


L  '  Y1_1. 


„  y  -  O'  -  l)v  _  'f  -  P  -  1 

(H'-v)*"1  (V  -  1  f  -1 


(9h) 


Solution  of  Eq.  (9)4)  for  Y  gives 

Y?-1  „  _ SU2  -  jp 


n(2-  •»  +  ^_2[i  +  (2->Vl  - 


[Y  -  (»-  l)v] 
(f  —  v)^  — 


(95) 


This  equation  gives  the  velocity-travel  relation  for  the  interval  after 
the  powder  is  burned.  . 


The  pressure  may  be  easily  found  as  a  function  of  the  velocity  and 
travel.  From  the  equation  of  motion,  Eq.  (8l),  and  Eqs. (86) ,  (89)  and 
(90),  it  follows  that 


12m'  dV  12  m  'Vu  dv  P.  dv 

P  =  -  V  —  =  - 2  v  —  =  _Q —  v  _ 

A  dX  AX,  dY  IX  dl 


4Xb 

Combination  of  Eq.  (92)  and  Eq..(96)  gives 


P  =  R 


r^- v) 

1  Y 


(96) 


(97) 


Thus  the  travel  and  pressure  corresponding  to  any  preassigned  velocity 
may  be  found  from  Eq.  ('95)  and  Eq.  (97)  ,  respectively.  Unfortunately, 
it  is  n't  possible  to' express  the  velocity  in  terms  of  the  travel  con¬ 
veniently  so  that  muzzle  velocities  must  be  found  by  successive  approx¬ 
imations  . 


FART  III.  CALCULATIONS  PROCEDURE 

1 14 .  Calculation  of  pressure-travel  and  velocity-travel  curves 

To  illustrate  the  method  developed  in  this  report,  we  calculate 
the  pressure-travel  and  velocity-travel  curves  for  a  hypothetical 
recoilless  gun. 


/ 


-  18  - 

The  gun  constants  are-  assumed  to  be  as  follows: 


D  Diameter  of  bore 

U. 13U  in.  (105  mm) 

A  Area  of  bore 

13.7  in? 

A^  Area  of  throat 

9.U3  in? 

vc  Volume  of  chamber 

336.9  in? 

Travel  to  muzzle 

80  in.  . 

PQ  Starting  pressure 

5000  lb/in? 

The  powder  constants  are: 

.  F  Impetus 

381*000  ft  lb/lb 

7  Ratio.- of  specific  heats 

1.223 

^  Covolume 

2 6.93  in?/lb 

p.  Density  of  solid  powder 

0.0396  lb /in? 

k  Nozzle  coefficient 

0.00398  sec-1 

B  Burning  constant 

0,00033  (in./sec)/(lb/in?) 

The  loading  constants  are: 

M  .  Weight  of  projectile 

29.76  lb 

W  Web  thickness -of  powder  grains 

0.031U  in. 

■  C  Weight  of  powder  charge 

.  8.379  lb 

First  calculate  the  following  preliminary  quantities: 

Quantity 

Numerical  Value 

Xo  -  vc/A 

26.03 

Ao  =  °Ac 

0.0231*8 

a  =  R  -  1  /p 

Pod  -  A0 JP)  ■ 

Jo  A0(l2F  +  aP0) 

10.17 

0.02769 

A 

c(bM 

o.U 

0.6303 

r  -  V(>  -  i)X  +  T 

1 .0680 

..  _  M  +  (1  -  ?l)C/3 

32.1  7I4. 

0.9370 

m'  =  I.OUm 

0.9933 
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Quantity  Numerical  Value 


r  =  (i  +  p)(?  -  i)  +  i 

1.3122 

„  (1  -  Jo)A 

b  m 1  (B/W) 

1250 

Sb  =  (7  -  Dm'Vb/2CF 

0.075U5 

g0  =  tO  ~  30)h 

0. 306U 

gi "  i  -  g0r 

0.6728 

□2  =  g0(7  -  l)/2 

0.03U16 

g3  =  i  -  -  gQr/r 

2 

0.5737 

Tb  /“g2  +  Vgg  +  g1g3\ 

T0”(  gx  ) 

0. 76)42 

+,  Tb 

>  .  1/2  r2  2  t0 

0.9058 

eavg'  1  -  S/f 

0.30)42 

Ji  *  i  -  wr 

0.2782 

_  R^avg  ~  "*//3 
a 

-0.8)438 

, ,  a  -  jpo  -rvyg) 
a  2lT  -  jQ) 

0.02)463 

CFm'  (B/\f)a 

el  - 75 - 

A 

1.953U 

CF(B/W) 

63  ~A 

699.1 

e3  =  I2m'e|/vc 

16.355 

r  =aA0j:[j'e1 

-0.1095 

^  =  Ao//*  +  a  Ao  jo 

0.)4006 

q  =  3o/ (cl)2  ei 

0.1832 

U  =  i(T  -  1 )  -  k^ex 

0.1080 

zb  =  vb/e2-  '  ‘  • 

1 . 7880 

Corresponding  to  any  value  of  the  parameter  Z  between  zero  and  Z-^ 
there  is  a  travel  given  by  the  equation* 


-  20  - 


where 


=  J(l  -  cl)  -  rS  +  ot  +  rZ. 
o 

The  pressure  and  velocity  for  this  travel  are  given  by 

p  _  ®3(n  +  Z  —  uZs) 

J(l  -  oO  -  rS 

and- 

V  =  e2Z. 


The  maximvsn  pressure  P^  occurs  at  Z  , 

ZP  -  f1  +§]/(1  +2“> 

or  Z,  ,  whichever  is  smaller.  It  will  be  observed  that  P  ,  the  maxi- 
b  p7 

mum  pressure,  is  used  to  obtain  a  value  of  Z  ..  However,  rP  /e3  is 

P  P 

small  in  comparison  to  unity,  and  thus  a  sufficiently  accurate  value 
of  P  can  be  found  by  a  method  of  successive  approximations  in  two  or 

ir 

three  steps. 

The  method  of  finding  ^bAQ  and  P^  will  now  be  given  in  detail. 

The  calculations  for  any  other  point  would  be  the  same  as  those  given 
below,  except  that  Z^  would  be  replaced  by  the  value  of  Z  at  the  chosen 
point. 


Quantity 


Numerical  Value 


-  (A  +  buq  ~  l)/2q  -  0.10 $9k 

a  (V'l  '  +  Ituq  +  l)/2q  3.561i3 

*  l/CriVT+TTuq)  9.0866 

=  l/(raVl  +  Uq)  0.17300 

-  0  -r^)^1  (1  +  r2Zb)~ta  U.U#7 

-  {j+(l-r1Zbrti-(l-r1Zb)[1+(l+r2Zbrta]}/[2r1(t1  +  l)3‘ 

h.609 


£/  J  and  S  are  the  functions  of  £,  u  and  Z  defined  in  Sec.  12. 
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« 


» 


xb 

Xo 

Pb 


L. 


’xb 


Quantity 

f 

-  J(l  -  «)  -  rS  +  <x  +  rZ^ 

_  e3(q  +  Zb  -  uZb)  . 

"  J(1  -  ex)  -  rS 


xb  \ 
■ - ' 


V 


Numerical  'Value 


3.3802 

83714 

62.00 

1  2.30 


The  cal  dilations  after  the  powder  is  all  burned  are  as  follows: 

Quantity  Numerical  Value 


Kb/C  =  gi  -  g0Mb/To 

1  =  2(gi  -  N^/C) 

N'/C 

Y  =1/^  +  1 


V?(2  -  fxw 


0.3223 

2.175 

1.U598 

0.9688. 


by 


The  travel  corresponding  to  any  velocity  V  greater  than  Vb  is  given 


t  s  7'  /y  523 

L  xor  r 

v  o 


where 


(7  -  i)v/vb] } "4i 
(V-v/vb)?_1  j 


The  corresponding  pressure  is  given  by 


P  »  R 


)(1  -  v/vb)^ 

I 


The  velocity  corresponding  to  any  travel  after  the  powder  is  all 
burned  cannot  be  found  directly.  The  method  used  here  is  to  calculate 
the  travel  at  several  velocities,  and  then  interpolate  to  find  the 
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velocity,  at  the  muzzle.  The  solution  to  our  example  corresponds  to 
taking  V  =  1355  ft/sec,  Y  =  1.209!?  and  L  =  80 .U  in.  The  correspond¬ 
ing  pressure  P  is  5010  lb/in? 

The  pressure-travel  and  velocity-travel  curves  for  this  example 
are  given  in  Fig.  2.  There  are  small  errors  in  these  curves  introduced 
by  the  mathemati  cal  approximations  explained  in  the.  text.  The  magni¬ 
tude  of  these  errors  is  shown  by  a  comparison  with  the  dotted  curve  in 
Fig. 2.  This  dotted  curve  was  obtained  by  a  point-to-point  numerical 
integration  of  the  fundamental  equations.  The  velocity-travel  curves 
are  indistinguishable  on  this  graph.  It  is  seen  that  the  mathematical 
errors  introduced  by  our  approximations  are  negligible  compared  to  the 
probable  physical  uncertainties. 

PART  IV.  CHARACTERISTICS  OF  RECOILLESS  GUNS 

The  only  existing  recoilless  gun  that  we  have  considered  is  the 
German  7«5-cm  L. G.  Its  characteristics  are  shown  in  Table  I  and  the 
pressure-travel  and  velocity-travel  curves  that  we  calculate  for  it 
are  shown  in  Fig.  3.  This  gun  operates  at  relatively  high  pressures 
and  therefore  requires  a  comparatively  short  travel. 

Table  I.  Characteristics  of  the  German  7»5~an  Leichtes  Geschutz. 


Chamber  volume,  v_ 

120  in? 

Density  of  loading, 

O.63  gm/cm3 

Ratio  A/A^ 

1.U5 

Maximttm  pressure, 

29000  lb/in? 

Muzzle  velocity,  Vm 

1210  ft /sec 

Weight  of  powder,  C 

2.75  lb  ' 

Approximate  web,  W 

0.02  in. 

Muzzle  pressure,  Pm 

9000  lb/in? 

Weight  of  gun  in  action 

320  lb 

In  order  to  examine  the  possible  characteristics  of  recoilless 
guns,  we  made  a  large  number  of  calculations  on  hypothetical  105-mm 
recoilless  guns  having  various  ballistic  parameters.  These  charac-  ' 
teristics  are  shown  in  Table  II.  The  following  conclusions  may  be  drawn. 


iet  hods :  of  tills; 
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FIG. 4.  PRESSURE -TRAVEL  AND  VELOCITY -TRAVEL  CURVES  FOR  A  HYPOTHETICAL 
RECOILLESS  GUN » THE  EFFECT  OF  VARYING  THE  WEB  W  IS  ILLUSTRATED. 
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(i)  The  muzzle  velocity  is  quite  insensitive  to  the  burning  rate 
or  web  of  the  powder;  consequently,  it  is  insensitive  to  the  maximum 
pressure.  It  is  also  quite  insensitive  to  the  starting  pressure. 

(ii)  The  maximum  pressure  is  very  sensitive  to  the  burning  rate 
or  the  web  of  the  powder;  as  a  consequence,  the  maximum  pressure  is 
very  sensitive  to  the  initial  temperature  of  the  powder.  This  short¬ 
coming  is  reminiscent  of  the  corresponding  problem  in  conventional 
rockets.' 

(iii)  It  Is  possible  to  design  a  recoilless  gun  operating  at  low 
pressures  that  has  a  very  flat  pressure-travel  curve.  This  is  not 
feasible  for  recoilless  guns  operating  at  high  pressures  since  the 
density  of  loading  cannot  be  increased  sufficiently.  This  is  shown 
in  Fig.  it. 

■  (iv)  The  muzzle  velocity  is  not  sensitive  to  the  cross-sectional 
area  of  the  throat  provided  the  web  of  the  powder  is  adjusted  so  as 
to  give  approximately  the  same  maximum  pressure. 

(v)  It  is  desirable  to  use  a  very  high  density  of  loading  in 
recoilless  guns. 

It  seems  probable  that  most  of  the  tactical  uses  of  rockets  can 
be  served  by  recoilless  guns.  The  weight  of  the  projector  and  a  few 
rounds  of  ammunition  is  quite  comparable  in  the  two  cases  since  the 
recoilless  feature  makes  possible  a  very  light  mount.  In  addition, 
the  recoilless  gun  has  the  definite  advantage  that  it  fires  convention¬ 
al  projectiles  with  good  exterior  ballistic  characteristics.  The  accu¬ 
racy  that  is  claimed  for  existing  recoilless  guns  is  quite  comparable 
to  that  for  conventional  howitzers. 

There  are  several  problems  that  must  be  met  in  the  design  of  a 
satisfactory  recoilless  gun.  First,  a  satisfactory  method  for  trapping 
the  powder  must  be  found  if  the  velocity  dispersion  is  to  be  kept  small. 
Second,  it  is  desirable  to  design  a  blow-out  disk  that  holds  to  a  fairly 
high  pressure  and  then  shatters  into  small  fragments  which  do  little 
damage.  Third,  it  is  -desirable  to  reduce  the  tremendous  blast  that  re¬ 
sults  from  the  escape  of  powder  gas  to  the  rear  of  the  gun.  The  blast 
in  the  Geraian  7*£~cm  gun  is  so  severe  that  the  personnel  must  use 
special  ear  protectors. 
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K 
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U) 

1  0  1 

o  1 

1  I 

H 

pb(AXb/Mi)r. 

Nozzle  coefficient;  k-  ? 

Atvf'‘/T0  (1  ~\/Tb/T”) 

2 

20(67(0  (1  -  nq/c) 

L 

in. 

Total  travel  of  the  projectile. 

M 

lb 

Weight  of  the  projectile . 

m 

slug 

H  +  (1  -  A)  C/3 

32 . 1 7U 

m' 

slug 

Effective  mass  of  the  projectile, 
ml  =  1 . 0l;  (M  +  9avg  C/S)/g  =  I.Olm. 

N 

lb 

Weight  of  powder  burned  at  any  time . 

N» 

lb 

Weight  of  powder  gas  remaining  in  the  gun 
at  any  time . 

Nb 

lb 

Value  of  N1  at  the  instant  the  powder  is  a! 
burned. 

No 

lb 

Powder  burned  before  projectile  starts  to 
m  ove . 

n 

mole/lb 

Number  of  moles  of  gas  formed  from  unit 
weight  of  powder. 

P 

lb/in? 

Average  pressure  in  the  gun  at  any  time. 

p 

xb 

lb/in? 

Pressure  at  the  instant  the  powder  is  all 
burned. 

po 

lb /in? 

Starting  pressure. 

ID 

xp 

lb/ in? 

Maximum  pressure. 

Px 

lb/ in? 

Pressure  on  the  base  of  the  projectile. 
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Tavg 

Tb 


T1 
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hi 

t2 

u 

V 

Vb 

v 

Yc 

W 

X 


ft  lb/ (mole  °K) 


°K 


K 


JK 


ft/ sec 
ft/ sec 


in? 


m. 

in. 


(N0/C)/ei(^)? 

Gas  constant  per  mole. 


a  A  j  j  e  . 

0dl'J2  1 


(\/l  +  Uuq  -  l)/2q  . 

(Vi  +  liuq  +  l)/2q  . 

Jy/dZ/j. 

Temperature  of  the  powder  gas  in  the 
chamber. 

Average  temperature  in  the  chamber  during 
the  gas  discharge.  Assumed  equal  to 

Kt;  ♦  ib). 

Value  of  T  at  the  instant  the  powder  is 
all  burned. 

Isochoric  flame  temperature  of  the  powder. 

T0/[>  -©(•?-  1)]. 


J 


1  + 


k  A+ 


(7  ~  D 


C  B/W 
1  /rxVl  +  liuq  . 

l/r3Vl  "+  Iraq  . 

|(7  -  1 )  -  ko01  . 

Velocity  of  the  projectile  at  any  time. 

Velocity  at  the  time  when  all  powder  is 
burned. 

v/V 

Volume  of  the  chamber,  equal  to  AXQ  .  • 

Web  thickness  of  the  powder  grains. 

Effective  distance  frcm  the  breech  to  the 
projectile  such  that  AX  is  the  volume 
behind  the  projectile. 


i 


;v 


j!  \  pj  f "  |  a  o  c  f:  f ™ F r> 


Effective  length  of  the. powder  .chamber*  equal 
to  vc/A. 

x/xb. 

(x/xj,-  a.. 


Value  of  Z  at  -which  maximum  pressure  occurs. 


(V>)  +  =>ao‘VS- 


lb/in? 


Ratio  of  heat  loss  to  kinetic  energy  of  projectile 
at  muzzle. 

/(/  -  1  )A  +  1  - 

Ratio  of  specific  heats 5  7  ~  (nR+  Cv)/CvI=1  +1 .9869  n/Cv. 

Psoudo  ratio  of  specific  heats*  J-  1  +  (l  +  /3)(>  -  l). 

Density  of  the  powder  gas  at  any  time* 

A  =  N 1  / [AX  -  (C  -  N)//°]  • 


lb/in?  Initial  density  of  loading*  equal  to  C/vc. 

-  Parameter  depending  on  C/k  and  usually  having  a  value 

a  little  larger  than  3*  see  A-lIj.2. 

F_+i?bF 

2T2  2T0. 

in?/lb  Covolume  in  the  equation  of  state  for  the  powder  gas. 


_  t 

C  (B/Yf) 

Average  value  of  S  during  the  burning* 


9  =  1  - 

avg 


C(B/W) 


Io/io  "t 


c  (  b/w-) 
kA. /C  (B/W). 


\J  e  t‘  w  is  6  fc..  <kJ 
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. _ e3Zb _ 

*  ~~  Pb^?(2  -?)(Xb/X0) 

p  lb/in?  Density  of  the  solid  powder. 
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ABSTRACT: 


A  system  of  interior  ballistics  is  developed  for  guns  in  which  the  recoil  is  eliminated  by 
the  rochet  action  of  the  powder  gas  flowing  through  a  venturi  in  the  breech.  Ballistic 
equations  are  derived  on  the  basis  of  simplifying  assumptions,  and  the  accuracy  of  the 
approximations  is  checked  by  comparing  a  solution  with  the  results  of  a  numerical 
integration  of  the  fundamental  equations.  The  system  is  Illustrated  and  applied  to  numerous 
examples  including  the  German  7.5-cm  Lelchtes  Geschuetz. 
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